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1. Introduction 

The motion of vortex sheets with surface tension has been analyzed in the set- 
ting of irrotational flows by Ambrose [T] and Ambrose & Masmoudi 12] in 2D, and 
by Ambrose & Masmoudi [3j in 3D. With irrotationality, the nonlinear Euler equa- 
tions reduce to the Laplace equation for the pressure function in the bulk, and 
the motion of the vortex sheet is decoupled from that of the fluid, thus allowing 
boundary integral methods to be employed. In a general flow with vorticity, the full 
two-phase Euler equations must be analyzed; in this situation, the vortex sheet is a 
surface of discontinuity representing the material interface between two incompress- 
ible inviscid fluids with densities p+ and p", respectively. The tangential velocity 
of the fluid suffers a jump discontinuity along the material interface, leading to the 
well-known Kelvin-Helmholtz or Rayleigh- Taylor instabilities when surface tension 
is neglected. The velocity of the vortex sheet is the normal component of the fluid 
velocity, whose continuity across the material interface T(t) is enforced. In addi- 
tion to incompressibility, the continuity of the normal component of velocity across 
T(t) is a fundamental difference between multi-D shock wave evolution, wherein the 
velocity of the surface of discontinuity is determined by the generalized Rankine- 
Hugoniot condition. Nevertheless, the problems are mathematically very similar, 
and we refer the reader to the book of Majda [6] for the analysis of multi-D shocks. 

In the incompressible, rotational flow-setting, very little analysis has been made 
of the two-phase Euler equations. With surface tension present, Shatah and Zeng 
[7] have obtained formal a priori estimates for smooth enough solutions, but the 
question of existence of smooth solutions remains open. In this paper, following 
the methodology of Coutand & ShkoUer [4] , we prove well-posedness for short-time 
for this problem. 

Let 17+ and fl~ denote two open bounded subsets of R'^ such that = 17+ U il~ 
denotes the total volume occupied by the two fluids, and F = il+ n fl~ denotes the 
material interface. We assume that it is the region fl^ that intersects dft. 

Let T] denote the Lagrangian flow map, satisfying 

?7t(a;, t) — u{r]{x, t),t) V a; G f2, t > , 
77(x, 0) = a; . 

Let n+{t), fl~{t) and r{t) denote r;(i)(f7+), r]{t){fl-) and r/(t)(r), respectively, and 
let andp^ denote the velocity field and pressure function, respectively, in Q^{t). 
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The incompressible Euler equations for the motion of two fluids can be written as 
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in 




(l.la) 


div 


= 


in 




(i.ib) 


[P]± 


= aH 


on 


r(t), 


(l.lc) 


[u ■ n]± 


= 


on 


m, 


(i.id) 


u~ ■ n 


= 


on 


on, 


(Lie) 


u(0) 


= Uo 


on 


{t = o}xn, 


(l.lf) 



where the material interface T{t) moves with speed u{t)^ ■ n{t), and p~ are the 
densities of the two fluids occupying and il^(t), respectively, H{t) is twice 

the mean curvature of r(i), cr > is the surface tension parameter, and n{t) denotes 
the outward-pointing unit normal on 917+ (i). 

Theorem 1.1 (Main result). Suppose that a > 0, T is of class i?^, Sfi is of class 
H^, andu^ € H^(^l^). Then, there exists T > 0, and a solution (u^ [t) , (t) , Q,^{t)) 
of (Hip with M± e L°°(0,r;ir3(n±(t)), p± e L°^{0,T;H'^-^{n'^(t)), and T{t) G 
i?^. The solution is unique ifu^ e H^-^{n^) and r e 

The paper is organized as follows. In Section [21 we establish the notation to 
be used throughout the paper. In Section [3] we establish low- regularity trace theo- 
rems of the normal and tangential components of vector fields with divergence 
and curl structure. In Section IH we introduce a regularized version of the Euler 
equations (jl.ip : the transport velocity and the domain are regularized using the 
tool of horizontal convolution by layers that we introduced in [4] . Additionally, a 
nonlinear parabolic regularization of the surface tension operator is made in the 
Laplace- Young boundary condition (|4.1t l) . Section [5] is devoted to the existence 
of solutions to ()4.ip . In Section [6l we obtain estimates for the velocity, pressure, 
and their time derivatives at time t = 0. Section [7] provides the pressure estimates 
that we need for a priori estimates. In Section [H we establish the K-independent 
estimates for the solutions of the K-problem (|4.ip ; this allows us to pass to the limit 
as the regularization parameter k — s- and prove existence of solutions to (jl.l|) . In 
Section [9l we provide the optimal regularity requirements on the data. Finally, in 
Section [10] we prove uniqueness of solutions. 

2. Notation 

Let n :== dim(f}) 2 or 3. We will use the notation H''{n+) (iJ*(f]")) to denote 
either i7"(17+;R) (iJ''(r2-; R)) for a scalar function or iJ''(17+;M") {H^n-;R")) 
for a vector valued function, and we denote the _ff'*(r2^)-norm by 

= and HwH^ _ = . 

The _ff*(r)- and iJ*(9r2)-norms are denoted by 

\w\s = \\w\\H-{r) and \w\s,dn = \\w\\H-{an) ■ 

For simplicity, we also use ||w||s,± and \uj\l .j. to denote + ||w^||s,- and 

\w^\l + \w~\l, respectively. That is, 

\w\l^^\w+\l + \w-\l 
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For s > 1.5, V^(T) denotes the space 

{w e L\0,T-,L^{n+)) w e L^iO,T;H'in+)nL°°iO,T;H'-^\n+))^ 



with associated norm 



kllvj(T) = sup ||w(i)lli/-i-5(n+) + / Ik (s)IIh=(o+) 

tG[0,T] 



where w can be either vector- valued or scalar- valued. The space Vi (T) is defined 
slightly differently, namely 

V1(T) = |w e L^{Q,T:L^{n-)) w E L'^{0,T; H^n-) n L°^'{0,T; H'-\n-))^ 

with norm 



\\w\\vUT}^ sup + / ||w(s)||^,(fj-)ds. 

te[o,T] 



As in [4], the energy function is defined as 

n 



2 

0,± 



+3.5,+ 



(2.1) 



3=0 



We use the notation = + + k to mean that 

J+=g+ + h+ , and r^9' + k-. 

3. Trace theorems 

The normal trace theorem which states that the existence of the normal trace of 
a velocity field w G relies on the regularity of divit (see, for example, [S]). If 

divw e H^{iiy, then w ■ N, the normal trace, exists in H~'^-^{dil) so that 



\\w ■ NWh-o.^oq) < C ||w||i2(f^) + II divwll^i(f^), 



(3.1) 



for some constant C independent of w. In addition to the normal trace theorem, 
we have the following. 

Theorem 3.1. Let w £ L^(ri) so that curlw e H^{Qy, and let ti, T2 be a basis 
of the vector field on dfl, i.e., any vector field u can be uniquely written as u"Ta. 
Then 



\\w ■ TaWn-o-^dn) < C ||w|||2(t^) + II curlwll^i(f^), , 
for some constant C independent of w. 
Proof. Given V' e H^ '^dQ), let 0^, G H^fl) be defined by 

A(j)a =0 in , 

= X Ta)tp on dQ . 

Then 



a = 1,2 



(3.2) 



(w ■ Ta)lljdS 



an 



j 

Jn 



cml w ■ (j>adx — / cm\(l)a-wdx 
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and hence 



Jan 



{wTa)ipdS < C ||u;||^2(Q) + II curlu;||Hi(o)' HaWmin) 

I L J 

< C ||w|| r2m> + II curlwll f/irnv || V''l|HO =(an) 



which imphes the desired inequahty. 

Combining (|3.ip and (|3.2p . we have the following: 



for some 



\\w\\H-o ^dn) < C*! l|w'llL2(a) 
constant C independent of w. 



II divw||^i(f^)/ + II curlw||ai(f^)/ 



□ 



(3.3) 



4. The regularized k-problem 



Let n' be an open subset of so that il+ CC fl' CC ft. In the following 
discussion, we will use M+ : H^-^{il^) H^-^{il) to denote a fixed bounded 
extension operator (from the plus region to the whole region) so that Af+w = in 
n"" for all V e i75-S(f7+). 

Let ti+ be the Lagrangian velocity in the plus region f2+, and Ve = M+w+ with 
Vk defined as the horizontal convolution by layers of Ve- Let ry„ = Id + u„(s)c?s be 
the Lagrangian coordinate (or fiow map) of v^,, and the Jacobian J7k, the cofactor 
matrices and the normal are defined accordingly. 

The smoothed K-problems is then defined as 



rje =ld+ / Ve{s)ds 
Jo 







in [0, T]xn^ , 

in [0, T]xn^ , 
in [0, T]xn^ , 
on [0, T] X r , 
on [0, T] X r , 
on [0, T]xdn, 
in {t^O}xn^ . 



(4.1a) 

(4.1b) 
(4.1c) 
(4.1d) 
(4.1e) 
(4.1f) 
(4.1g) 



where = ?7k,i ■ f7K,2 is the induced metric on F. Note that since M"*" extends v'^ 
continuously to the whole domain Q, rj'^ ^ Vk vt ^ Ve on F. 

Remark 1. Since AI+v = in fl'" for all v e H^-^{n), = Id and n,^ ^ N on 
dQ. Therefore, the boundary condition J^. JJ fJ can also be written as ■ N ~ 
where N denote the outward pointing unit normal of Q~ on dQ. 



5. Existence of solutions for the regularized k-problem 

5.1. The iteration between the solution in ri+ and il^. Let {v~^,v~q) E 
V+^(r) X Vl-^(r) X Vi^'^(r)/R be given, and let v+ be the horizontal convolution 
by layers of v~^. Define Ve = A/+u+, the extension of with the associated 
Lagrangian map = Id + Jq Ve{s)ds and cofactor matrix 0,^ = i7k(V^k)~^ where 
Sk = det(V7yK) is the Jacobian. The normal vector is then defined by 
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The process of finding solutions to (|4.ip consists of finding solutions to the fol- 
lowing two problems. First, in the plus region fJ"*", we solve 

P^J.wl + (a,)^rj = in [0, T] x r!+ , 

(a«)>^=0 in [0,T]xr!+, 

r = 9- (7X3(77) • n„ - kAo(w • n„) on [0,T]xr, 

w{Q)=u+ on {t = Q}xVt, 

where w = m+ o 77^, r = p+ o 77^, 77 = Id + w(s)ds and Ag = ^a[^/gogo^^ dp]. 
Once the solution (?«, r) to (|5.ip is obtained, then in the minus region f2~, we solve 



(5.1a) 
(5.1b) 
(5.1c) 
(5.1d) 



Jnvl + {d^YAv 





+ Miq^j 


= 


in 


[0, T]xVt- , 


(5.2a) 






= 


in 


[0, T]xVl- , 


(5.2b) 




V ■ 




on 


[0,T] xr, 


(5.2c) 




V ■ 


= 


on 


[0, T]xdn, 


(5. 2d) 




v{0) 


= "0 


on 


{t = Q}xVL, 


(5.2c) 



where 7; = w~ o 77^, q = o fj^^. 

This process introduces the map $ : («+, v^,q) t-^ {w, v, q), and the fixed-points 
of <I> provides solutions to problem (|4.ip . 



5.2. Estimates for the solution in 17+. The only difference between (|5.ip and 
the one phase problem studied in [4j is the presence of the term q in the boundary 
condition (|5.1b ). We note that if q is smooth, then by exactly the same argument 
as in [4], the solution to ()5.ip will be also be smooth, depending on the regularity of 
the initial velocity Uq . Therefore, for q given in L^{0,T; H^-^{V,~)/M.), we replace 
mt) by 

r = q^- aLg{ff) ■ n^, - kAq{w ■ ti^) on T , (5.3) 

where q^ denotes the horizontal convolution by layers of q. The solution w"^ and r*^ 
to (|5.1b .). (|5.1b ). (|5.3p and (jS.lti l are smooth functions satisfying 
rt ^ 

11.5^+ + ^^Iw" ■ n^\l^ ds<N{ua) + C{K,v+,q)^ft, (5.4) 



where C{k,v~^ ,q) denotes a constant depending on k, v^-°'(T) i II9|Iv^^(T)- Note 
that although this constant depends on p'^ as well, we omit this dependence in the 
estimate since it is a constant. 

The divergence and curl estimates as in 4J can also be carried on so that 

\\cuvlw^i^^+ + \\divw'\\l,,+ < Niuo) + C{k,v+) f \\w^h^+ds (5.5) 

Jo 

for some constant C(k, v^) independent of the smooth parameter e. Estimates ([57 
and (|5.4p imply that 

ft , rt 



t 



\\w'{s)\\i^^+ds < -N{uo) + C{k,v+) 



and the Gronwall inequality implies that 



"'0 



\w^[s')\\l^^^ds'ds 



\\w\s)\\l^^+ds<C{ti,u+,v+,q)Vt. 



(5.6) 
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By studying the elliptic problem for r'^ with the Dirichlet boundary condition 
we find that 

\V{s)\\l^^+ds<C{K,u+,v+,q)^t. 
57)) implies that wl £ L^{0,T; H'^-^{n+)) and by interpolations, 



sup \\w'ml+ < \\u+\\l+ + CiK,u+,v+,q)VT. 
te[o,T] 



further implies that 



||r^||2,+ < N{uo) + C(«, 6, u+,v+,q)Vi + S\\q\\ 



(5.7) 
(5.8) 
(5.9) 



It also follows from (15.1b ) that \\wf\ 



y2.5(-2.-) shares the same bound as ||?''^||y3.5(-y-), 



I.e., 



|w^?llv^-=(T) 



< N{uo) + C(k, S, u+, v+,q)Vi + S\\ 



9llt 



(5.10) 



These e independent estimates enable us to pass e ^ and obtained solution {w, r) 
to problem (|5.ip with estimate 



lkllv=-(T) + ll«^tllv-=(T) + Ikllv3-=(T) < Niuo) + C.,sVT + S\\q\\i^ , (5.11) 
where Ck,^s is the short hand notation for C{k, S, u'^ , ,q). 

5.3. Estimates for the solution in fl^ . We will set up a iterative scheme in 
order to show the existence of a solution to problem (|5.2p . Let = J~^{aKYy 
For a given w e Vl-^(T) with wt G V^-^(T), we solve first 

^f[^kjifc=-P"^'<.^?^^^:fc-P~4'[^'(^^~^ -^^e"-''K,],fc in (5.12a) 

= ^P^N • + - • + - ^^fr^)'^!£^K] on r, (5.12b) 



= p [w- n^t 



on dfl. (5.12c) 



Once a solution to (|5.12p is obtained, use this solution q in (|5.2b .) and solve the 
transport equation 



P 



in Q 
in 



W(0) = 

Suppose we can prove that v is actually in the space we start with, then a fixed-point 
of the map : w ^ v provides a solution to problem (j5.2p . 

We note that in this iterative scheme A is always fixed with estimates 

A{t)\Ur,,_<C{v+)Vt (5.13) 

for some constant C depending on ||w+||y5.5|.j,j but independent of k. Therefore, by 

assuming that T is small enough (so that C(w+)T is small), it follows from elliptic 
theory (see |S]) that 



\1\h.5.- 



< C 



l^llk- + lht|lk+ + ll^ll2.5,H 



\q\\ho.-<c 



W\ 



2.5 



- + \\wt\\u.+ + 



1.5,+ 



Combining these two estimates and (|5.1ip . by interpolations we find that 
lkllv3.5(T) < N{uo) + C^,sVt + S\\q\\l_ + Ct\\\w\\1,_,^^^ + ||iDj2..5(j,) 



(5.15) 
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For the regularity of v, we mimic the divergence and curl estimates as in [3]. In 

n-, 



(5.16) 



where 



a function of Vu, Vu and Vr/, where we use the identity A^A^^ ,, = AjA^ ,,. Let ( 
be the solution to 

i.e., C, is the flow map of the velocity field A'^{v — Ve), then 



curl Wo + / o Cds 



(5.17) 



Since 



([5T7)) implies 



if(C(y,s),s)ds oC\x,t)= / if(C(a;,s-t),s)ds, 



(5.18) 



We use (|5.18l) as the fundamental equality to proceed to vorticity estimates in il 



Since WCmlz.. <Mq + CT\\ 



C{v~), ((ET5)) implies that 



\\cnT^^v\\l^^_<C{v') N{u^)+ / \\v\\l^^^ds 
Transforming back to the domain /^^(^^ ), we find that 



(5.19) 



II curlu||j:^3.5(^^(n-)) < C(w ) N{uo) + \\u\\H3■^f^^(n-))ds 

Jo 

We remark here that the restriction of obtaining higher regularity is mainly due 
to the presence of VA in B{v) that comes from the transport term. Boundary 
conditions (\5.'2b ) and (15. 21 1) imply 

h • = \\u ■ ^||?f4(,^^(r)) + \\u ■ ^||l/4(f^^(af^)) < C||w||i5_+ . 

These two estimates and the divergence free constraint div u = lead to 



X{T) <C\\w\\^,^j,^+C{v-)\TN{uo)+ / X{t)dt 

J 

where X{T) ^ \\ 

^ll?f*^(ri„(o^))'^^' Therefore, the Gronwall inequality implies 

^0 







\Hm ^n.{n~))ds < [1 + Civ-)T]Niuo) + C^jVT + S\\q\\i 



or equivalently. 



Ikllis -dt < [1 + Civ-)T]Niuo) + C^^sVt + d\\q\\ 
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For T even smaller (so that C{v )T is small), it follows from (|5.2b .) that 



\\vt\\i^^_dt<C ||t'|lk- + ll9ll3.5, 



ds 



< 7V(mo) + C^^s^T + 5\\q\\i_ + CT 



klivl-S(T) + \\wt\\\,2.^T) 



Therefore, 



||w||y4.5(y) + ||l't||v2.5(-y) + ||9||v3-5(T) 



< N{uo) + C^,s^+5\\q\\l_ + CT \\w\\l.j.^r^) + WwtWl^..^^) 



(5.20) 



In the following sections, we will always assume that the initial input q satisfies 
||g||^3.5^j,^ < N{uq) + 1. We can choose a fixed but positive S < ^^^^ — )"-i-T) ' ^^'^ 
let L be the collection of those elements v E L^{0, T; H'^-^{il^)) so that 

||w||v4.5(j.) + ||wt|lv^-5(T) < ^(^o) + 1 • 

For a fixed k > 0, we choose T small enough so that 



C^.sVf+CT 



N{uo) 



1 



1 

< - 
- 2 



Clearly the map maps from L into L. Similar to the proof in Section r5.41 5* can be 
shown to be weakly continuous in L^(0, T^; H^-^{n~)). Since L is a closed convex 
set in L^(0, T; i/^-^(f7~))), by Tychonoff fixed-point theorem, there is a fixed-point 
V of the map which provides a solution to (|5.2p . Uniqueness follows from the fact 
that (15.211 is linear. 



(5.21) 



Remark 2. It follows from H5.20\) that 

Ik llvl-5(T) + Ikt|lv2-5(T) + ll9llv3S(T) < ^("o) + 1 



5.4. Weak continuity of the map Let {v^, qm) converges weakly to (w^, q) in 
the space L^{Q, T; i75-^(r!±))xL2(0, T; H^-^{n-)/R), ^v+, y-, qm) = iw„,,Vrn,q^) 
and $(«"'", g) = {w, v, q). Suppose that J7^, a^m^ ^tm constructed from 
accordingly. By the property of convolution by layers and the weak convergence, we 
have (i/„„, o^m' '^Km 



) converges to {J^,a^,n^) strongly in [L°°{0,T; H^-^{n^))f 



Since {wm,rm) satisfies 



0+ 



0+ 



qn 



I \ In 



and {wm,r,n) are uniformly bounded in L^{0,T; H^'-^{n+)) x L^{0,T; H^-^{n+)), 
it follows that there exists {w,f) so that 



q + Lg{rt) ■ nt\ {dt)lN,^^dS G i/i(f^+) . 

By the uniqueness of the solution to the linearized problem, w = w. Similar 
argument shows that the solution {vm,,qm) to problem (j5.2p with all the fixed 
coefficients constructed from converges weakly to {v,q), the solution to problem 
(j5.2p . Therefore, the weak continuity of the map <i> is established. 
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5.5. The fixed-point argument. The only thing we need to check is that if there 
is T > and a closed convex set K C L^{0^ T; H'^-^{n+)) xL^{0, T; H^-^{n-)/R) x 
L'^{0,T; H^-^{n~)) so that $ maps from K into K. Let K be defined as the 
collection of those elements {w^,q) G L^{0,T; H'^-^{n+)) x L^{0,T; H* '^{n-)) x 
L^{0,T;H^-^{n-)/R) so that 

ll^^llv=-=(T) + WwtWv^^^iT) < + 1 , 

lk"|lv4.5(T) + l|wr|lv2.5(T) + h\\v^'^(T} < ^(""o) + 1 • 

Recall that 6 is fixed from the previous section. Similar to the proof in the previous 
section, we choose T > small enough so that 

C^jVT + T{N{un) + l) < \ . 

Then by estimates (|5.1ip and (|5.2ip . the map $ indeed maps from K into K. 
Therefore, the Tychonoff fixed-point theorem implies the existence of a fixed-point 
(w, q) of $. 

Remark 3. This T is k- dependent. 

Remark 4. Once a solution to problem (^TTp is obtained, without loss of generality, 
we may assume that the pressure and its time derivatives satisfy the Poincare in- 

Jn- ' 

determined up to the addition of a constant (constant in space), we can replace 
and by q^ — g(= Q^) and q^ — q{= Q^) 

||g||g,± = iiQii^ < ciivgii^ = c\\vQ\\l^ . (5.22) 

5.6. Estimates of the divergence and curl of the velocity field. 

5.6.1. Divergence and curl estimates. In f2^, we can apply exactly the same tech- 
nique as in i4j to conclude the following lemma. 

Lemma 5.1 (Divergence and curl estimates in Let Li = curl and L2 ~ div, 

and let 770 :— ri{0) and 

:= P{\\u^\\2.5+n,+ , |r|3.5+n, Vk||Uo ||l.5+3n,+ , %A:|r|i+3n) 

denote a polynomial function of its arguments. Then for j — 1,2, 

n+l 

sup ||V^ij?7+(t)||2 5+„_+ + ^ sup ||ijC»tS+WII?.5+n-fcH 



te[o,T] ' ■ ^te[o,T] ''^ 

n+l T 

^ / II V^i,a^+||L+„-fe,+di < A/o+ + CTV{ sup E^{t)) , 



(5.23) 



fe=0 



Similar to the way of obtaining (|5.19p . the following lemma is valid as well. 

Lemma 5.2 (Divergence and curl estimates in Let n, Li and L2 be defined as 

those in Lemma\5A\ and 



^(11^*0 l|2.5+n ,-, |r|3.5+n, V^\\Uq ||l.5+3n \/K|r|i+3n) . 

Then for j = 1, 2, 

n+2 

J2 sup \\L,d^v-{t)\\l,+,,_^^_ < + CTV{ sup E,{t)). (5.24) 
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5.6.2. H '^■^ -estimates for on the boundary T and dQ. By (|4.1b ). 
(curl 1-4)' = e^jk [[S^j - {anYj]v^uj " (aKttfjV^t " (a^tYj^^t'^t 



div4, ^{6^- A'^vfl, ~ {Aut)lvf - i{Att)lvt:^ - i{At) 



f7± 



Since vf^^ G i (0, T; _ff (with k, dependent estimate), cnv\v^^ and divu^ 

are both in L'^{Vl^) and hence by p.3p . j:/-o.5(aQ±) exists. For (^9 G 



dn- 



dn- 

Since = Id and — Q outside 17', we find that 

Jan- 



— £ijk 

= 0, 



dn 



where we use the boundary condition (|4.1b ) and (|4.1f ) with = on F to 
conclude the last equality. Therefore, 



curl • ifdx 



< C 



t\vttt\-o.5,±+V{E^{t)) M 



which implies 



cnr\vttt\\HHn-Y < C t\vttt\_o.5.± + V{E,{t)) 



Similarly, 



div V. 



< C 



t\\vttt\\H-o-^(an±)+V{E,{t)) 



Therefore, by ([O]) . 



< C 



k*'ttllL2(ji) + II divt)j=^J^i(n)' + II curl 



ttt\\m{ny 



<CT\\v^tt\\H-o-^an±) + CV{E,{t)). 
It then follows from choosing T > small enough that 

\vut\-o.5,± + \vut\-o.5,an < CV{E^{t)) . 



(5.25) 
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6. Estimates for velocity, pressure, and their time derivatives at 

TIME t = 

In this section, we estimate the time derivatives of the velocity and pressure at 
the initial time t = 0. We use Wk, k = 1,2, 3, and qt, £ = 0, 1, 2, to denote dfv{0) 
and 9(^(0). Let ip^. be defined by 

[J-^ai{J'^a'lip^),k],j ^0 in n+ , (6.1a) 

(f^ = -aLg{r]e) ■ - kAo{v ■ Uk) on T , (6.1b) 
(fK. — Q on dil , (6.1c) 

and the quantities ipQ, ipi and 1^92 be defined by (Pk{0), ^Kt{0) and (p^^ttiO), respec- 
tively. 

Let Qq and Qq denote the initial pressure q{0) in 0+ and f2~, respectively, then 
and q^ satisfy 

-^A(g+-(^o) = /+ in n+ , (6.2a) 

-^Aq-=r in n-, (6.2b) 

^^ = -<-iV on dn+, (6.2c) 

;^^ = ("^^+^(-(o)-'c7)^o)-^^ on on-, (6.2d) 

where = (Vm^)^ : Vuj^, and denotes the unit normal of F from into 51+ , 
or the outward unit normal of dfl. 

Remark 5. The right-hand side of i6.S!b ) is in fact f^ — (v^ (0) — Uq) ■ Vdivwp 
while the last term is zero by the divergence free constraint of the initial data. 

For all tp e H^{n+) n H^ifl^) so that -0+ = on F, we have 

^ / \7{q+ - ipo) -"^^dx + — ( \/qf^ ■\/'i/jdx = [ f+ipdx + [ f~^/jdx 
P Jn+ P Jn- Jn+ Jn- 

{w^ ■N~S/^-,^.-u^)^dS- [ {w+ +V^o)- N^dS . (6.3) 

Since [v ■ n^]± = and w+ • = on dfl, it follows that 

w+ ■ N + u+ ■ n^t{0) ^ ■ N + Uq ■ n^tiO) on F, 
• = on dn, 

and hence (|6.3p implies 

[ \/q^ ■ \7i>dx + — I V(q(7 - (po) ■ S/iljdx = / f+tpdx + [ f^^jjdx 
P Jn+ P Jn- Jn+ Jn- 

+ f[iu+ -u^)-n^t{0)+Vipo- N]HS+ f V^-,^.-u--NHS. (6.4) 
Jr Jon " 

Let Qo ~ q^ — ipo in and Qq = q^^ in . Since Q'^ = on F, we can 
use Qo (and its difference quotients) as a test function in (|6.4[) . Since n„t(0) = 
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—Qq^ {uqk,^P ■ N)ld_a and ||we(0)||fc__ < C||uo||fc,+ , by standard difference quotient 
technique, for s > 1.5 for n = 2 or s > 1.75 for n = 3, 

< CV{\\uo\\l±) + C(|r|2^i.5 + k\uo ■ N\l^,,,) . 

By the ehiptic estimate for ip in (|6.ip together with (|4.1b ). we find that for s > 
l + s(n). 



l«^ill'-i.± + ho\\l± < cr{\\u„\\l^, |r|^+i.5, \\^u+\\ 



S+2.+ ) 



(6.5) 



For j = 1,2, the quantities qf^ and satisfy 



1 



+ 2{j - l)Vgo±(v4,V4, - VwfjN + Br 



in n"^, 

on dfl"^ 



where 



- 4(Vu^Vwi«) : (Vuonf ~ 2(Vw^Vmok) : (V?«i«)^ + i^w^^f : Vm^ 
+ div[2(Vu;i«)^Vq± + 2(VMo«VMo«)^Vg± - (Vu;i«,)^Vg±] , 

fcl = [ - lig,,(«e(0)-^' - V)Vfc + («e7(0) - «^r)^0,5- + (^^e'''(0) " "^)u;i":^- 
^2 = [(24kj"L/ - «^Lj)(«e"-''(0) - + («e«(0) ^ 

+ 204,(0)- VX 



and 



01 = -2Vuok : V^<^o - V<^o • Awok , 

02 = -4(ViioKViioK) : V^(y3o + 2VwiK : V^<^o - div(VMoKVuoK) • V(^o 

+ A?i;i« • V^o + 2Auok ■ V^i + 4Vmok : V^i + 2V[(Vwok)^V(Pi] : Viio« , 
Si = {y-{^) - UofWu„^{Wu^fN ~ {v;M ~ wYfWuoN ~ {v-{0) ~ UofWw^N, 



Bo 



Vii(, (2VuokVmok - Vwik){v^ (0) - Wq ) + Vmq (Wett(O) - ) 
+ Vw2"(ve"(0) - ) - 2(Vwo Vuo«)(w~j(0) - w^) - 2Vw^Vu^^{v-{Q) - u^) 
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where uqk, "m^Ik and W2k are Af+w,!^, M^w^^ and M+w^^, respectively. Similar to 
the estimate of , since on F, 

q^-'fj=qj forj = i,2, 

[i;+(0) - vrAO)] ■ N = -2K]± • n,t{0) - [uo]± ■ fi.tt(O) , 
[v+M - v^M] ■ N ^ ~3[W2]± ■ n,t(0) - 3[wi]± • n^tt{0) - [uo]± ■ n^ttt{Q) , 
and on dfl, 



we find that for s > 2, 



('2||s-2,± + \\H\\1-1,± 



< CV{\\uo\\U,\T\i^,,, 
and for s > 3, 

Ikalls-s.i + h2\\l-2,± 
< CVi\\uo\\l±, |r|2+i.5 + k\uo ■ N\l^,,, + K\wt ■ Af|^+o.5 + «k2+ • A^ls-o.s) : 
where we also use the boundedness of the extension operator M so that 

||9,^e(0)|U.-<C||u,+ ||^,+ . 

7. Pressure estimates 

The estimates for the pressure and its time derivatives are exactly the same as 
(12.1) in [4]. In [4], the L^-estimate for the pressure is found by studying a Dirichlet 
problem, but in the two-phase problem with fixed outer boundary, the L^-estimate 
is not necessary because of the Poincare inequality. Therefore, 

lkWll3.5,± + htmL,± + kttmU < cv{E,{t)) (7.1) 

for some constant C independent of n. 

Remark 6. The estimates for , and q^^ require the control of ||u^||3.5,_, 
ll't'r Il2.5,-; ll^^tt lli-5.- '^'^'^ ll^tttllo,-j respectively. This is the only reason we need to 
include the estimates for d(v~' into our definition of energy Note that we do 

not need ||77~|j4.5__ in order to control dfq^ . 

8. k-independent estimates 

We also make use of the following inequality which follows from Morrey's in- 
equality (see (2.6) in [4]). For U € ^^^'^'(r), 

\U^{x)-U{x)\<CK^''i\DU\p (8.1) 

Test (HT|) against a function e iji(f7+) n i7i(f7") with (/j" • = on d^, 

n+ Jn- L ' J 

ia,)iq+ip+'dx + f {a^)lq-^ip-'dx = . (8.2) 
n+ ' Jn- 

Similar to those estimates in the K-independent estimate consists of studying 
the three time differentiated problem, three tangential space differentiated problem, 
and the intermediate problems with mixing time and tangential space derivatives. 
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Most of the estimates are essentially the same as those in [4 , and in the following 
sections we only list those terms which required further study. 

Before proceeding, we remark that those energy estimates in [4] can be refined 
a bit further. For example, the energy estimate for the third time-differentiated 
K-problem ((12.6) in [4j) can be refined as 



sup 

te[o,T] 



WvtttWl + \vu ■ n^\l + / \^d^v n^Wdt < Mo{5) + 5 sup E^{t) 

te[o,T] 



+ CTV{ sup E^{t)) + C{5) \\vt\\l^ + \\v\\l,^ + hellis + / \\Vlivu\\l^dt 

tG[0,T] L 

where the difference is not having 



C sup 

te[o,T] 



P{\\vA\l.,) + P{Ml, + P(hllL) + cp{\\^v^A\l 



on the right-hand side of the inequality. To see this, for example, one such term 
comes from estimating 



sup \P{v,dr],,)\L^{dn) j \^/K^'^v■n^\l\^fK^fv^,\2dt. 

tG[0,T] Jo 

Since P{v,drii^)t £ L°°{0,T; L^{T)), by the fundamental theorem of calculus, 
sup \Piv,d7j^)\L^rr)<Mo + CTV{ sup E^{t)) 



te[o,T] 



te[o,T] 



and hence 



JQ 



sup |P( 

te[o,T] 

<d sup E^{t) + CTV{ sup E^{t)) + C{S) [ W^v+Wl^dt, 
te[o,T] te[o,T] Jo 

instead of having || VK'^tt|li2(o t-h2.5(q)) in the bound shown in 4J. Therefore, the 

energy estimates we cite from :4 will have only one polynomial type of term in the 

bound: CTrisnpt^[o,T]E^it)). 

In this section, we will make use of the following equality which follows from 



■ {vttt - -Vut) = n^ttt ■ {v - v+) - Sn^tt ■ {vf - ) - 3n«,t • («+ (8.3) 



8.1. Estimates for the third time-differentiated K-problem. Three time dif- 
ferentiate (18. 2|) . and then use vttt as a test function and integrate in time from to 
T, we find that 



E 

sign— ± 



Jn^'s" 



Vflf^^dxdt 



lo Jn- 



(a«)-(i' '-v,')v^ 



Vffldxdt = . 
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1=1 P 

Jn- 

i-T f 



Jo Jn+ 



v^f-ldxdt 



ttt 
■T 



v^fldxdt . 



ttt 



lo Jn- 



ttt 



The worst terms of Xi is when all the time derivatives hit v^k, while the other 
combinations are bounded by CV{E,^). Therefore, 



Xi < 



Jo Jn 



p-{a.)](v-= - v-=)v-l,vuldxdt + CTV{E^) 

p-{a,)){v-^ - v-^)\vttt\\dxdt + CTV{E^) 

, ^ p-(a,)^^ («-■'■ - v-')Nk\vttt?dSdt + CTV{E, 
~ Jo Jon- 

The boundary of 11^ consists of F and dVL. On dVi, — Id and Ue 
by 6If), 



JO Jn- 



0. Therefore, 



Jdn 



p-{a^f^{v-^ -v-^)Nk\vut?dS ^ f p- {y- ■ N)\vtu\^ dS ^ . 

Jdn 



_ 1 

On r, since = f+ and — gn ^{aK.)'jNk, boundary condition (|4.1b ) implies 
that 



^p-ia^^{v-i - v;i)Nk\vttt\^dS ^ J^^p-[v ■ n^]±\vttt\^dS = 0. 



Therefore, 



Ii < CTV{E^) 



(8.4) 



The worst terms of T2 is when all the time derivatives hit q. Therefore, 



l2< / {a^)iq+t.jVmdxdt 
Jo Jn+ 



^ ^ JO 

sign— ± 



Jn 

g^{a^iN,v^^dS 



{a.)iqut,j^mdxdt + CTV{E,) 



qtlrMlvlll'^dx]dt + CTV{E,) 



= I21 + 222 +CTV{E,). 
For X21, it follows that 



121 = 



Jo Jt 



+ +i _ - -i 

qttt'^ttt qttt'^ttt 



dSdt 



/g^Qtttivttt - Vttt) ■ n^dS + y/g^iq'^ - q )ttt{vttt ' 'n^,)dS dt . 
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By (|8.3p and substituting —aLgijie) ■ — kAo(w • ni^)nK for ((7+ — q ), we apply 
the estimates as in 14 to obtain 



I21 ^~ 1^ y/g^^lttt 



dSdt 



+ 5 sup E,{t)+M^{5) + CTV{ sup E^{t)) + C{5) \\v+\\l^^^ 



te[o,T] 
2 



tG[0,T] 



+ IK'+ll3.5,+ + hc=||l5,+ 



Integrating by parts in time, since 
the same techniques as in 4^ , we find that 



fg:{vt -vi)- n^tt )] ^ e (0, T; i2 (D) , 



3 



< 5 



JO Jv ^ ^ 

sup E,{t) + A/o(<5) + CTV{ sup E^{t)) . 
6[0,T] *e[o,T] 



9it 



gf,{v^ - ) • "-Ktt dSdt- / q^^^/g^{vt - ) ■ n^ttdS 
J t Jr 



t6 



using 

t=T 

t=o 
(8.5) 



Let the first and the second term of X2i„ be denoted by l2i„ and ^21^ i respectively. 
Integrating by parts in time, 



i21. 



Jo Jr 



+ — T' 



dSdt 



t=T 

I V^Qu^Kt ■ (vtt - Vu)dS 
Jr *=o 



By [Vg^n^t ■ iv+ - v^t)] ^ e L^O, T- H-^ %T)). Since q^, € L°°(0, T; H°-^{T)), 

it follows that 



^ Vg^QttnKt ■ {vtt - Vtt)dS 



* ^<(5 sup E^{t) + Mo{S) + CTV{ sup E^it)) . 
t=o te[o,T] te[o,T] 

Again by (|5.25p . we can estimate the first integral of 221^2 and obtain 

<S sup E^it) + Mo{S) + CTVi sup E,{t)) . (8.6) 

tG[0,T] tG[0,T] 

For l2iaj^ , integrating by parts in time again. 



^21„, = 



^ / <ltt V9^{v^ - V )■ dfuf, + [y/g^{v^ - V )]tn,,ttt 
Jo Jr L 

r t=T 
- / y/g^qttiv^ ~v^)-nKtttdS 



dSdt 



t=T 
t=0 



The second term of X2i„j can be bounded by CTP (supjgjQ E^{t)) since the inte- 
grand is in i°°(0, T\ L^{T)). Since n^ttt ~ Fi{drif^)dvKtt + F2{drif^, dvf,)dv^t, by the 
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factthat ^^{v+ -v-){F^+F2{dr]^,dv^)dv^t) e L°°{0,T; L\r)) s^nd H°-^{T) 
^-0.5 (-p-j (duality pairing, 

< Mo + 



t=T 
t=0 



Mo + CTVi sup E^{t)) |g,7(T)|o.5 |af««(T)|_o.5 + 1 
te[o,T] J L 

< Mo + Mo||q„(T)||i,_[||«„tt(T)||i.+ + 11 +CT7'( sup E^it)) 

-I te[o,T] 

< Mo{S) + S sup E,{t) + CTV{ sup , 

t6[0,T] te[o,T] 

where ||u„tt(T)|lo,+ < Mq + CT7^(supjg[o,T] -^k(^)) ^-^id Young's inequality are used 
to obtain the last inequality. 

It remains to estimate the first term of X2iaj^ in order to complete the estimate 
of l2i- We write the first term as 

- r I {q+ - qu)Vr.{v^ - v^) ■ dfn^dSdt X3) 
Jo Jr 

T 







By n^t = -g'^'^ivn.a ■ ni^)rii^,j3, 



QttVd^i'"^ - V ) ■ 111^,09"'^ {vKttt.a ■ n^)dSdt 



JT 
T 







where the second integral is bounded by CTV {sup f^^o,T] E^it))- For the first term, 
QuVd^iv^ - v^) ■ VK„0gZ^{vK.ut,a ■ n^)dSdt 



dSdt. 





It follows from H°-^{T)-H~^-^{V) duality pairing and (I5.25P that the term with 
{v^ttt ■ n^^a) is also bounded by CrP(suptg[g -j.] E^it)). 

Let ^ be a non-negative cut-off function so that supp^ C IJ^ suppa^ and ^ = 1 
on r. Integrating by parts in space, since dT = (f), by the divergence theorem, 

q:^t^/gZ{v^ - v^) ■ ■qn^pg'^^^Vt^tu ■ n^)^adSdt 





T 



Jr 

T 





T 



-1' )-f7K,/35« qtt{vnttfn^,)dSdt {=Iii) 



In 



dxdt , 



7) 



/o Jn+ L 

where w is an i/^-^-extension of to il^. By (|5.25p . X41 can be bounded by 
Cr7^(suptg[o T] ^K.{t)) as well. For the rest terms, there are two worst cases: when 
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the derivative dj hits ^ or w^ttt- For the latter case, by inequaUty 

WiMKttt.j - [^{a^)l'"ltt,M\G,+ < CK\\a^\\3^+\\vtu\\i,+ ■ 
This inequahty together with the "divergence free" constraint impHes 

mcin)HtuJo.+ < C4a,h,+ \\vut\\i,+ + Cr{E,{t)) 
and therefore by Young's inequahty, 

lo Jn+ 

<5 sup E,{t) + CiS)TV{ sup E,{t)). 
te[o,T] te[o,T] 

For the former case, we make use of the equation ()4.1b ) to substitute {aK)iq,k for 
vl- Therefore, in this case the worst term is 



Jn+ Jo Jn+ 

Let Qi — {a^)lq^^ ■ and = i.iv'^ — w) ■ rji^^^g"^ . By the definition of horizontal 
convolution by layers, we find that 

/ / dc.Q,Q,^F''dxdt = y^ I I {^aQ^){ei)[p*hP*h{Q^{Oi))]F"{ee)dydt. 
Jo Jn+ ^ Jo -'[0,1]=' 

Since {dc^Q.m ^ eZd^iQM), 

i-T I- 

[^o.Q^){el)[p*hP*h {Q0t))]F-^{ei)dxdt 



"'[0,1] 



1 



2 Jo J[osY 



^^\p*hiQ^m)\'F''{ei)iee)ldydt+ / Rdt, 



where R — p^h 



F''i9e)iQi)l^^Q^i9i) - F'^{et){Qt)lp^h d^Q^iOt) 



and by in- 



equality dHU), since VQ, - Fi{d-q^)vut + ^2(i9?7„, dv^,dv^t)^q + F^idr]^, dv^)Vqt, 

[ \R\dt<Cnl ||i^(0)e,||H.i,oc([o.i]3)||a(Q,(0))|U2([o,i]3)dt 
Jo Jo 

<Mq{5)+5 sup E,{t) + CTV{ sup E^{t)) . 

te[0,T] t6[0,T] 

Integrating by parts in space, 

/ ^^\p^H{Q^{6^WF^{el){Q,)ldydt<CTV{snv E^{t)). 
lo J[o.iY te[o,T] 

Combining all the estimates above, we find that 

Ii<MQ{5) + 5 sup E^{t) + C{5)TV{ sup E^{t)) . (8.8) 
te[o,T] te[o,T] 

Now we turn our attention to X22 before estimating X3. By the "divergence free" 
constraint (|4.1b ). 



I22dt= 

sign— ± 



JVP'sn 



(n Vii^'Si" I 0/ y sigm o/„ V„signi 



dxdt . 
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As shown in P], it follows from integrating by parts in time that 

/ / {a^tt)iv^;qutdxds < S snp E^{t) + Mo{S) + CTV{ snp E^{t)) . (8.9) 
JQ Jn± t6[0,T] te[o,T] 

For the first and the third term, we follow Q and obtain 



E 

sign— ± 

^ E 

sign— ± 



sign 
Ittt 



JSl^isn 

T 



dxdt 



JSl^isn 



s sigm O S sign* 



E 

sign— ± 



Jn='s> 



s signs O S sign* 



qllTdxdt (=122 J 
qllTdxdt (=l22j 



+ 5 sup i?«(i) + Mo(<5) + CT7'( sup £;,(i)) + C(<5) + + ll^+lls 5 + 

te[o,T] te[o,T] L 

+ he||4.5,+ + / WV^vttWh. + dt 

Jo 

Using the "divergence free" constraint again, 



sign— ± 



JOsiBn 



<^ sup E^{t) + Mo{5) + CTVi sup E,{t)) , 

te[0,T] t6[0,T] 



(8.10) 



where we apply estimates similar to (|8.9p again from [3]. 

Integrating by parts in time (and space if there is VKttt or Vttt), since a„ = Id on 
dfl and Vk, — outside fl' (or near d^l), we find that 



2'22. < 



Jr 



Jr 



N.dSdt (=J22,J 
NjdSdt {=122 J 



+ S sup E^{t) + Mo{S) + CTV{ sup S,(t)) , 
te[o,T] te[o,T] 

where similar estimates for the lower order terms are obtained as those in [3]. It 
follows from ([g?^ and that 



121,, < / Km\^o.5V{E^{t))dt<CTV{ sup E^{t)) 
Jo «e[o,T] 



122.. = -3 



Jr 



fg^J^, ^{a^Ys (vm - «ttt) • n^^Qu + (9+ - q )tt{vut ' n^) 



< S sup E^{t) + Mo{S) + CTV{ sup E^{t)) , 
te[o,T] te[o,T] 



(8.11) 
(8.12) 



where we use the boundary condition (|4.1b ) in the second term and apply the same 
estimates as in [4]. 
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For X3 , we use the boundary condition (|4.1d ) in X3 and obtain 



/o Jr 

= ^31 + ^32 



aAg{rie) ■ + kAo{v~^ ■ Uk) \/g^iv^ - v ) ■ dfn^dSdt 



The worst term of X3 is when the time derivatives hit the highest order term. Since 

T 



"tttlll.5,+ 



'^tt II2. 5,+ 



I32 < CTV{ sup E^{t)) - 
te[o,T] 

Integrating by parts in time, 







dt < Ek{T), by Young's inequality, 

'5|l%/^4JL,+ + c('^)llV^4ll2.5,+ 



dt. (8.13) 



131 = 



Jr 
a 

r 



Ag{rie)-ni^ Jj;:{v^-v ) ■ d^n^dSdt (^Jsi^) 
. ittt 



(^^3lJ. 



For X3i^ , it foUows from integration by parts (in space) that 



< 



cr.g'^*(Wt^ 7) • nny/g^iv^ ~V ) ■ ?7K,/3ff"'^ (wKtt,a<5 ' n^^jdSdt 

Jt 

+ CTV{ sup E^{t)). 
te[o,T] 

By the definition of the inequahty above imphes that 

l3U < - [ I <y[d^iP*h vtt) ■ n,]F''-<\ds{p^h vtt) ■ n^UdSdt 
Jo Jr 

+ CTV{ sup E,{t)), 
te[o,T] 

where F°''^^ = y/gZg"^ g^^ {v^ ~ v^) ■ rj^^p. Since F"')''^ is symmetry in 7 and (5, it 
follows from integration by parts that 







cr{d^vtt ■ n^)FZ'\dsv+ ■ n^)dSdt + CTV{ sup E,{t)) 

te[o,T] 



and hence 



(8.14) 



2'3i„ < CTV{ sup E^{t)). 

te[o,T] 

Integrating by parts in space, the worst term of T^n^ is 

Since F"''^ £ L^(0, T; i°°(r)), integrating by parts in space, we find that 

l3U<S snp E^{t) + C{S)\\vt\\h^+ + CTr{snp E^{t)). (8.15) 



te[o,T] 



te[a,T] 
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Combining all the estimates above, 



httt\\l,± + \vtt ■ n\i + iV^Vut ■ n-Klldt 



sup 

te[Q,T] 

<Mo{6)+5 sup E^{t) + Ci6)TVi sup E^t)) + C{S)\\\v+\\i^^+ (8.16) 

t6[0,T] *6[0,T] L 

+ Ik+ll3.5,+ + he||1.5,+ + / Il\/^4ll2.5, + rf^ 



We also need controls for \vf^ ■ n\i. It follows from inequality (|8.ip and the 
fundamental theorem of calculus that 

\w-{n-n^)\i<CK\Mo + CTV{ sup |«;|i|?7+|4.5 • 

L te[o,T] J 

Therefore, by (|4.1b l and the fundamental theorem of calculus, 

Wt ■ n\i < \v^t ■ "kIi + Wt ■ {n - n^)\i 

< IKt - • '^kIi + • ^^kIi + ■ (" - 

< |2(wt+ - Uj") • n„t + (w+ - w") • n^ttli + l^tt -"-11 + l«tt ' (»^ - »^k)|i,± 

< Mo{S) + S sup E^t) + CTV{ sup + |4 • n|i . 

te[o,T] te[o,T] 

Having this additional inequality, we find that 



kttt|lo,± + k« • "li.± +/ \^vtu ■ n^,\ldt 



sup 

te[o,T] 

<Afo(^)+<5 sup E,{t) + C{6)TV{ sup + [||«+||i5^+ 



te[o,T] 



k+ll3.5.+ + lh+ll4.5.+ 



te[o,T] 



|V«W«ll2.5,+'^* 



(8.17) 



8.2. Estimates for the second time-differentiated K-problem. Similar to 
(12.33) in [4J, let £,ddl act on (|4.1b ) and test against £,dvtt, we find that for Si > 0, 



sup \d^Vfn\l^+ \V^d\+ ■ n^\l ±dt < Mo{Si) 
te[o.T] Jo 



(8.18) 



+ Si sup E,{t)+C{di)Tri sup E,{t))+Cidi) / \\^v+\\l^_^dt. 

t6[0,T] t6[0,T] Jo 

8.3. Estimates for the time-differentiated K-problem. Let £,d^dt act on (|4.1b ) 

and test against ^9^wt, we find that for 62 > 0, 



sup \d^v-n\l± + [ \V^d\+ ■ n^\l ±dt} < Mo{62) (8.19) 
te[o,T] 'Jo ' ^ 

+ S2 sup E^{t)+C{52)TV{ sup E^{t))+C{52) [ \\V^v+ \\l ^ _^dt . 
t<E\0,T] telo.T] Jo 
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8.4. The third tangential space differentiated K-problem. Similar to (12.37) 
in [3], the study of the boundary condition (|4.1t i') leads to the following important 
elliptic estimate: 

sup |V^77+(t)l5± < Mo + C sup E,{t) + CTV{ sup E,{t)) . (8.20) 
te[o,T] ' te[o,T] t6[0,T] 

Let £,d^ act on (|4.ip and test against £,d^v, by (|8.20p . we find that for ^3 > 0, 

sup |5V • ^lo ± + / iV^d^V^ ■ n^\l ±dt 
te[o,T] Jo 

<Mo{S3)+S3 sup E^{t)+C{S3)TV{ sup E^{t)) . (8.21) 
tG[0,T] te[o,T] 

8.5. A polynomial-type inequality for the energy and the existence of 
solutions. Combining the div-curl estimates (|5.23p . (|5.24p . the energy estimates 
([8T7| . ([8A81) . ((8J9| . ([OT]) . we find that 

E^{t) < AM6,di,62,S3) + iS + 6iC{S)+d2C{di)+ 630(62)) sup E^^t) 

te[o,T] 

+ C{6,6i,62,63)TV{ sup S«(t)). 

te[o,T] 

Choose 6 > and 5^ > small enough so that 6 + 6iC{6) + 620(61) + 630(62) < ^, 
then the inequality above implies 

sup E^{t) < Ma + OTV{ sup E^{t)) . (8.22) 
te[o,T] te[o,T] 

Therefore, there exists Ti > independent of k so that 

sup E^{t) < 2Mo. (8.23) 

tG[0,Ti] 

This K-independent estimate guarantees the existence of a solution to problem p.ip 
by passing k — * 0. 

8.6. Removing the additional regularity assumptions on the initial data. 

In the previous sections, we in fact assume that v is smooth enough so that we can 
directly differentiate the Euler equation (|4.1b ) and test with suitable test functions. 
This requires higher regularity of the initial data, namely, S H^'^-^{fl^) and 
r e iJ^. As in [41, this can be achieved by mollifying the interface by the horizontal 
convolution by layers and mollifying the initial velocity by the usual Fredrich's 
mollifiers. 

8.7. A posteriori elliptic estimates. As in by exactly the same proof, we 
find that for T sufficiently small, 

sup [\T{t)\5.5 + \\v\U.5,± + \\vth,±] < Mo , (8.24) 
te[o,T] 

where Mq is some polynomial of AIq. 
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9. Optimal regularity for the initial data 

In the previous discussion, the existence of the solution requires the initial data 
G H'^-^{n^). We show that this requirement can be loosened to £ H^{fl^) 
and r e i/^ in this section, by assuming that we already have a solution to the 
problem. 

In this section, we study the problem in the Eulerian framework. To start the 
argument, we define the energy function £{t) first. Let £{t) be define by 

^(t) = |r(i)l4 + ll"^llH3(f2+(t)) + \\u ||H3(o-(t)) + 11"*^ llffi-5(n+(t)) 

+ ll"t llHi =(0-(t)) + \Wtt\\L^(n+(t)) + ll"ttllL2(0-(t)) • 

Then for the pressure function p^, we have the following estimate: 

\\P^\\H^-Hn+{t)) + lb"ll?/2 =(o-(t)) + \\pt\\m{n+{t)) + \\Pt\\m{n-(ty) < CV{£it)). 

The estimates for curlu^ are essentially identical, while the estimates for divu^ 
are trivial because of the divergence free constraint (|f .lb ). Therefore, 



sup 

te[o,T] 



< Moid) 



curlw+||^2.5(o+(t)) + II curlw ||H2.5(o-(t)) + || curlw+||^i(n+(i)) 

div li" 



curliii ||^i(o-(t))ll divw+||^2.5(o+(t)) 



divWf 



div 



\m{n-{t)} 



(9.1) 



6 sup £{t) + CTV{ sup £{t)) 
te[o,T] te[o,T] 



where Mo{S) = Mo{\r\l \\u+\\l^, \\u^\\l_,6). 

Remark 7. The reason for not analyzing the problem in the ALE formulation is 
that in the minus region, the transported velocity {v^ ~ only as regular as 

Vrf^ , which is less regular than the velocity . This prevents from obtaining the 
estimates for curl?;^ in H'^-^{n^). With the Eulerian formulation, the transport 
velocity is H^{n^{t)), and the analysis goes through. 



Note that the need of that 77 is more regular than u (or v) is only for the study 
of the K-problem, where the estimate of the boundary integrals with artificial vis- 
cosity K requires that 77^ is as regular as ^/kv (which is "more regular" than v by 
the definition of the energy function Ek). This observation implies that without 
worrying about the K-terms, the energy estimates still follow. Therefore, by the 
identities 

f fiy,t)dy 

and 



/ (/t + V„±/)(2/,t)dy 
Jn±(t) 



f{y,t)dSy^ / fir]{x,t),t)y/gdSx , 
r(t) Jr 

we can show, as shown in the previous sections, that 



sup 

te[o,T] 



IKWhin+it)) + ll"ttlli2(o-(t)) + |5 " • "lo,± + \9vt ■ n|o,: 



< Mo{S) + S sup £{t) + CTV[ sup £[t)) , 
te[o,T] te[o,T] 



(9.2) 
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where the interior estimates are for the Eulerian velocity while the boundary 
estimates are for the ALE velocity v^. 



In addition to |r(i)||, it suffices to establish bounds for |9 



'2„± 



lffo-5(r(t)) 



and 



\duf^ ■ '7i|^2(p(()), where m denotes the unit outward normal of 17+ (i). Wc remark 
here that we use different notations to distinguish the "normal" on F and the normal 
on T{t). In general, n = m o rj. 

The bounds for \duf^ ■ "^|^2(p((')) follows from the energy estimate (|9.2p . Since 



±2 i 



on r(t) 



multiplying — (^J^f^^ ° both side, by \\5 — a||2,+ ^ C'(t) we find that 

\do.uf ■ m|i2(r(t)) < C\\d^v^ ■ n\o + \dJv^'v^/y\o\ + CTV{ sup £{t)) 



te[o,T] 



< Mo ((5) + 5 sup £{t) + CTV{ sup £{t)) . 
te[o,T] te[o,T] 



(9.3) 



For the bound of • w|^o.5(r(t))' estimate • n|Q 5. Similar to 

the a posteriori estimate in [3], by studying the boundary condition 



where rj and g are formed from , we find that 

\d\+ ■ n\l < MoiS) + S\v+\l + C[v{\r\l„£mrj+ - Id\l + \p+ prio 
\d\+ ■ n\l < Mo{S) + S\v+\l + C[v{\r\l„£mrj+ - ld\l + \pt-p7\l 

and hence by interpolations, 

sup \d'^v+ ■ n\l ^ < MoiS) + S sup £{t) + CTri sup £{t)) 
te[o,T] te[o,T] t6[0,T] 



C 



\\pt\\'H^(n+(t)) + \\Pt llffi(o-(t)) 



By the elliptic problem 



in n^{t), 



t_ 

dn 



(<+V„±u± + V„±uf)-n on r(i) 



we find that 

„±ll2 ^ r-' I1V7„,± . /'v7„,±\T|i2 I |L,± 



\\Pt\\m{n±{t)) <C W'^ut : (Vw ) |lffo.5(n±(t)) + + V„,u + V„wt||i2(j^±(t)) 

(9.4) 



<Moi6) + S sup £{t) + CTV{ sup f (i)) , 
te[o,T] te[o,T] 



where we use 



terms. Therefore, 



to estimate \\u 



ttllL2(0±(t)) 



and Young's inequality for the other 



sup \d'^v+ ■ n\l^ < Mo{6) + 6 sup £{t) + CTV{ sup £{t)) . 
te[o,T] te[o,T] te[o,T] 

By similar argument of obtaining (|9.3p , we find that 

sup \d'^u+ ■m\]jo.^n+(t))<Ma{5) + 5 sup £(i) + Cr7'( sup £{t)) . (9.5) 
te[o,T] «e[o,T] te[o,T] 

The estimate of \d'^u~ ■ '7i|^o.5(r(t)) follows from the boundary condition (jLlti 'l. as 
discussed in the previous sections. 
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It then foUows from dm]), ((121), and that 



sup 

te[o,T] 



^(i)-|rWlI <Mo(<5) + (5 sup £{t) + CTV{ sup . (9.6) 

t6[0,T] te[o,T] 



With this estimate in mind, we can estimate ||?'^llff2.5(Q±(i)') in the way we obtain 
(|9.4p and find that IIp*II/^2.5(-q±((')') satisfies the same inequality. Let h be the height 
function of T{t) over F. By exactly the same argument as in [4], 



sup |/i(t)|^4 < Mo((5) + (5 sup £{t) + CTr{ sup £{t)) , 
te[o,T] t6[0,T] te[o,T] 



and hence 



sup \r{t)\l < MoiS) + 6 sup E{t)+CTV{ sup . (9.7) 

tG[0,T] te[o,T] te[o,T] 

Combining (|9.6p and (|9.7p . by choosing (5 > small enough, we obtain the same 
polynomial-type inequality as (|8.22p . and therefore there exists a T > so that 

sup £{t) < 2Mo. 
te[o,T] 

This proves the claim of the optimal regularity of the initial data to obtain the 
solution to (jl.ip . 

Remark 8. The argument in this section can also be used to prove the existence 
theorem for the one phase problem studied in [A\, provided the same regularity of 
the initial velocity uq and the initial interface F are given. 

10. Uniqueness of solutions 

Suppose that {v-^,q^) and (w^, g^) are both solutions to (|4.ip (with k ~ 0, J = 1) 
with initial data Uq g H^{fl^) and F G . Let rjl and rj^ be defined as in Section 
|4](with associated cofactor matrices a^ and a^), and set 

m = \\m\\7 + E \\dtv'mi-i.,K± + E ii^^v wiii5-i.5.,± ■ 

fe=0 k=0 

By the existence theorem, both £i{t) and E-zit) are bounded by a constant Aio 
depending on the data uq and F on a time interval < t < T for T small enough. 

Let w = — v'^, We = M+w^ with associate flow map Q = Jq M~^w~^ds, and 
r = — q^. The goal in this section is to show that w = by showing that the 
energy function 

i?(0-ll«(i)ll3,± + ll«*WllU± + IK'«Wllo,± 
is actually zero for a short time. 

10.1. The divergence and curl estimates. In 0+, and satisfy 

p'^v:^^ + a'^jqj = for (v.q) ^ {v^ , q^) or (v^q) ^ {v'^ , q^) . 

Let EijkajVr act on both sides of the equality above and form the difference of the 
two equations, after integrating in time from to t, we find that 



p+ cnrlw+\t) = / (a' - a^)i{vl)\^ + [(a^)! - 5i]wtl 



ds 



e^Jk[[a%{t) - 5i]wy{t) + e,,, / [a' - a%{vl)y " {dtay.w^' 



ds . 
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Therefore, by \\a^{t) - <5(i)||3.5,+ < CT, 



sup II curlw+(t)|||_^ < CT sup ||w+(t)||3.+ 



te[o.T] 



te[o,T] 



where C depends on Aio only. By the "divergence free" constramt a^w j = 0, we 
similarly have 

sup ||divw;+(t)||2+ <CT sup ||u;+(t)||2 + . 



te[o,T] 



te[o,T] 



For the divergence and curl estimates in , let -D^ and 5^ denote the Lagrangian 
velocity in 17", that is, 

y = dtT]-' — o fj-' in , 
where is the Eulerian velocity in fJ^. The same argument as above shows that 



sup 

te[o,T] 



||curlw(i)||^_ + ||divw(t)||^_ < CT sup II w(t)||3 _ , (10.1) 

te[o,T] 



where w = — v'^. We now convert (jlO.ip to the inequality with w replacing w. 
Let C^' = (t/J)"! o (77^)-i and V = VC^' for j = 1, 2. Then 



II curlu;(t)|||+ = ^ \\e.,,[i' o - o C%\\i_ = ^ \\s^Mib')lv]^ " {b')li 

^\\l+ + \\m\\i 



1=1 



< CT sup 
te[o,T] 



where we use ||(^^ — ^^)(i)||2,- < C'^sup(g[o,T] 
time derivative of 6^ — b^. Similar argument shows that 



by studying the 



diviullg _ < CT sup 

te[o,T] 



3_<CT sup 

te[o,T] 



w\ 



3,- 



3,+ 



llilli 



Thus for r > small enough, we find that 



sup 

te[o,T] 



||curliy(<)||^_ + ||diviy(t)||^_ < CT sup ||w(i)||^_. 

te[o,T] 



The estimates for the divergence and curl of Wt are similar, so we omit here. In a 
nut shell. 



sup 

te[o,T] 



curlw(t)|||± + II divu-(<)|||± + II curlu;t(t)||2^5 ± 



divwt(i)||g5± <CT sup E{t) . 

te[o,T] 



(10.2) 



Remark 9. We cannot obtain estimate U0.2\) by studying the equations for v" 
(with the transport velocity) directly, since it also requires the study of C,^ as we did 
in the estimates of the K-problem. This requires e H^{n^) at least. 
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10.2. The boundary estimates. First we 


note that {w 


7") satisfies 








in 


[n Tl X D,"^ 


no Sa"! 


r l\j ±i 12 l\j 2±i 

(a = (a - a 




in 


[0, T\xn , 


(lO.dbj 




B 


on 


[o,r]xr, 


(10.3c) 


■ ni — ■ ni + b 




on 


[o,r]xr, 


(10.3d) 


W~ -711=0 




on 


[0, T] X 917 , 


(10.3e) 


w;=^(0) = 




in 


= 0} X 1]± , 


(10.3f) 



where 



2± 



b- 



(t;^+ - ) • ("-2 - ni) . 

The main difference between (|10.3p and the uniqueness argument for the one 
phase problem (see Section 15 in [4]) is on the additional term b. In order to obtain 
estimate similar to ()8.22p (except that in the uniqueness proof, we only study the 
second time differentiated problem), we need to estimate the integral 

i-T f 

rui^tt -Wtt)- nidSdt. 



By mM, 



tt - Wtt) ■n = btt- 2{w+ -Wt)-nt-{w+-w )ntt ■ 



The only term we need to worry about is the integral with integrand r^"^ (w^+ — 
v'^~) ■ (n2 — 'ni)tt- The worst term of this integral is (after integration by parts in 
time) 



^"^Kt•nl)<,-5^"^(^.^^i-2)<, 



dS 



+ 



la/3/ 1+ N 1 2a/3/ 2+ „\„2 



dSdt 

t=0 

dSdt. 

t 



Add and subtract terms to form the integrand in terms of ~ Ve^ "-i ~ '^2 or 

9^ — 9^ ■ By Young's inequality, the first term (time boundary term) is bounded by 
{5 + C{5)T) supig[Q E{t), where C{5) depends on 5 and A^o- For the second term 
(time interior term), the worse term occurs when time differentiating dvt- For this 
worst case, we can transform the surface integral to the interior interior using the 
divergence theorem as we did in (|8.7p . Therefore, 

g'"^{vlX ■ ni)illp - 9"^\vll ■ n^)^!.^ dSdt < {S + C{S)T) sup E{t) , 

■ J* ie[0,T] 

The estimates with the addition of the forcing F, the right-hand side of (|10.3b ). 
and B is already done in [3]. It suffices to show that \d'^w ■ nijo.s has the same 
bound. However, since 

\Bt\l^<CT sup E{t)+C\w\l,<{S + CiS)T) sup E{t) , 



te[o,T] 



te[o,T] 
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by study the first time derivative of (|10.3b ). similar to the a posteriori estimate, we 
find that 

sup \d'^wni\l^<{5 + C{5)T) sup E{t) . 
te[o,T] «e[o,T] 

Therefore, with (|10.2p wc conclude that E{t) satisfies 

sup E{t) <{5 + C{5)T) sup E{t) , 
te[o,T] te[o,T] 

which implies for T small enough, E{t) = and hence w — Q. In other words, we 
establish the uniqueness of the solution to the problem. 
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